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Abstract 

u 

Simple  approximations  are  given  for  the  average  run  lengths  of  CUSUM  tests:  both 
with  and  without  the  fast  initial  response  feature,  and  stuting  from  a  quasi-stationary  state. 
Numerical  examples  illustrate  the  accuracy  of  the  results. 

'■  ' 

1.  Introduction. 

For  1/  =  0, 1, 2,  •  •  • ,  =  oo  let  denote  the  probability  measure  under  which  yx, yj,  •  •  • ,  Vi/ 
are  independently  distributed  with  probability  density  function  go  and  are  indepen* 

dently  distributed  according  to  y,  which  is  assumed  to  belong  to  some  family  of  probability 
density  functions.  Let  gi  denote  one  particular  possibility  from  the  set  of  densities  g,  let 

=  log[ji(yk)/yo(yfc)]i  and  Sn  =  *1  ^ - hin-  A.  CUSUM  test  for  detecting  the  change-point 

u  is  defined  by  a  stopping  rule  of  the  form 

r  =  inf{n  :  S„  -  min  Sk>b}. 

k<n 

See  vanDobben  de  Bruyn  (1968)  for  a  background  discussion  of  CUSUM  tests. 

The  stopping  rule  r  is  traditionally  evaluated  in  terms  of  its  average  run  lengths:  J5oo(r), 
the  in  control  average  run  length,  which  should  be  large,  and  £3o(r),  the  out  of  control  average 
run  length,  which  should  be  small.  Poliak  and  Siegmund  (1985)  suggest  that  one  also  consider 
E„(r  - 1/ 1  r  >  */)  for  large  */,  i.e.  the  average  run  length  starting  from  a  quasi-stationary  state. 

Assuming  that  go  and  g  can  be  imbedded  in  an  exponential  family  of  distributions,  Sieg¬ 
mund  (1985,  Theorem  10.16)  has  given  approximations  for  Eoo{r)  and  £b(r).  The  main  result 
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of  this  note  is  an  approximation  for  Ei,(r  —  1/  \  r  >  v)  which  is  valid  when  1/  and  6  are  large. 
We  also  improve  Siegmund’s  (1985)  approximations. 

To  provide  a  “fast  initial  response”  to  detect  changes  occurring  at  1/  =  0,  Lucas  and 
Crosier  (1982)  modify  the  definition  of  r  by  putting  5o  =  i,  S„  =  *  +  ii  +  •  •  •  +  i„,  and 

T  =  inf{n  :  5„  -  min  St  >  fc), 

Oo^lc<n 

where  cq  =  inf{n  :  5n  <  0).  We  also  give  approximations  for  the  average  run  lengths  of  this 
modification  and  show  by  a  numerical  example  how  our  approximations  can  provide  guidelines 
for  selecting  the  value  of  x. 

A  more  precise  description  of  the  approximations  is  given  in  Section  2.  Section  3  contains 
numerical  examples  to  illustrate  the  accuracy  of  the  approximations. 

Poliak  and  Siegmund  (1985)  in  the  context  of  Brownian  motion  studied  the  Shiryayev- 
Roberts  procedure  as  a  possible  competitor  to  CUSUM  tests.  It  would  be  desirable  to  have 
approximations  similar  to  those  of  the  present  paper  for  that  process.  Poliak’s  (1987)  approx* 
imation  for  the  Poo  average  run  length  is  adequate,  but  otherwise  it  seems  a  difficult  problem 
to  develop  simple,  comparably  good  approximations  for  the  Shiryayev-Roberts  process. 

2.  Approximations 

It  will  be  convenient  to  give  our  results  in  a  canonical  form,  which  usually  is  not  the 
most  suitable  form  for  any  particular  application,  except  perhaps  for  tests  for  a  change  in  a 
normal  mean.  We  develop  here  our  notation  and  definitions  ab  initio.  They  are  not  necessarily 
consistent  with  Section  1. 

For  $  in  some  real  interval  containing  0,  let  F$  denote  a  probability  distribution  of  the 

form 

dP«(x)  =  exp(ffi  -  xlf{$)]dFo{x). 

Then  ^'(tf)  =  /  xdF$  and  V’"(^)  =  /  -  [^’(®)]*  >0-  By  an  affine  transformation  if 

necessary  one  may  assume  without  loss  of  generality  that  V’'(0)  =  0  V’**(0)  =  !•  Then 


»gn{9)  =  agnyp'{$). 


For  eiich  <  0  (^i  >  0)  assume  there  exists  >  0  (^J  <  0)  such  that  V'(^t)  =  (»  =  0.  !)• 

By  the  convexity  of  V*  the  value  of  ^5  (^i)  is  unique.  Let  A,-  =  Oi  -  6*^.  Note  that  Aq  <  0  <  Ai- 

Let  ^0  <  0  <  and  v  =  0, !,•••, +oo.  Let  Py  denote  the  probability  under  which 
are  independent  with  distribution  and  Xy+i ,  *  *  *  are  independent  with  distribution 
Ft^.  Of  course  Py  depends  also  on  Oq  and  $i,  but  we  suppress  this  dependence  to  simplify  the 
notation.  When  i/  =  0  or  oo  so  the  entire  sequence  is  independent  and  identically 

distributed,  it  will  be  convenient  to  write  to  emphasize  dependence  on  9.  (Note  that 
according  to  this  convention  Pq  now  denotes  Pt  for  9  =  0,  not  Py  for  v  =  0.) 

Let  0  <  X  <  b,  Sn  =  X  +  xi  -{ - 1-  and 

r  =  inf{:  —  min  S*  >  b), 

ao<k<n 

where  cq  =  inf{n  :  n  >  0, 5^  <  0}.  When  x  >  0  we  shall  write  Pf  to  denote  dependence  of 
probabilities  on  the  value  of  i. 

It  will  be  convenient  to  use  the  following  additional  notation.  Let 

r+(r_)  =  inf{n  :  Sn  >  0(5„  <  0)}. 

Put  Pi  =  .Eo(S'r^)/[2£b(5,+)]  and  po  =  ^{S’_)/[2£b(5,_)].  In  general  po  and  pi  must  be 
computed  numerically.  For  an  algorithm  see  Siegmund  (1985,  Theorem  10.55  and  Problem 
10.7).  In  the  case  of  a  standard  normal  distribution  pi  =  —po  ^  .583.  Also  let  yi  = 
jE^(5®^)/[3f^(5T+)]  -  pj  and  72  =  •Eo(S?_)/[3£t)(5T_)]  -  Pq-  Like  po  and  pi  the  constants 
7o  and  71  must  in  general  be  computed  numerically.  However,  their  role  is  comparatively 
minor  and  one  might  choose  to  neglect  them.  For  continuous  and  symmetric  Fo>  hence  in  the 
normal  case,  it  is  easy  to  use  the  Wiener-Hopf  factorization  of  the  distributions  of  Sr^  and  5r_ 
(e.g.  Siegmund,  1985,  Theorem  8.41)  to  show  that  71  =  70  =  E!o{x\)/l2. 

We  shall  assume  throughout  that  Fo  >s  strongly  non-arithmetic  in  the  sense  that 

limsupl  f  e*^*dFo(i)|  <  1. 

lAHoo  J 

Similar  results  hold  if  Fo  is  arithmetic  provided  the  definitions  of  po,  pi„  70,  and  71  are 
appropriately  modified. 
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Our  first  result  is  a  refinement  of  Theorems  10.13  and  10.16  of  Siegmund  (1985). 


Proposition  l.Let  $q  <  0  <  Bi  and  A,-  =  ff,-  —  as  above.  Suppose  0i  —*  0  and  b  —*  oo  in  such 
a  way  that  A, -6  converges  to  a  finite,  non-zero  limit.  Then 


(2)  Bg.(r)  =  [A,V''(ffO]~Mexpl-A<(6  +  fii  -  Po)]  +  A<(6  +  pi  -  po)  -  1}  +  7i  -  Tfo  +  o(l), 


Eo(r)  =  (b  +  Pi  -  Po)*  +  7i  “  7o  +  o(l)- 


Uniformly  for  b  bounded  away  from  0  and  1 


^<.(^)  =  lA.V>'(tf,)]  *{exp[-A,(5-l-  Pi  -  Po)]  -  exp[-A,(i  -  po)] 


+  A, (6  -I-  Pi  -  *)}  +  7i  +  o(l)» 


?=  (fc  -  *  +  Pi)(^  +  *  -I-  Pi  -  2po)  +  71  +  o(l). 


Proof.  The  expansion  (3)  is  obtained  by  Siegmund  (1975);  and  (2)  to  one  degree  less  precision 
is  given  by  Siegmund  (1985,  Theorem  10.16).  The  slightly  refined  form  of  (2)  given  here  may 
be  proved  by  following  the  approach  of  Siegmund  (1985)  in  conjunction  with  the  arguments  of 
Theorems  1-3  of  Siegmund  (1979).  The  details  are  omitted. 


To  obtain  (4)  and  (5),  note  that  for  all  x  and  0 


E^(r)  =  Ei{N)  +  Pt{SN  <  0}£^(r), 


where  N  =  inf{n  ;  Sn  ^  [0,6)}.  According  to  Siegmund  (1979,  Theorem  3) 


pnss  <  o> = + „(»-■), 

exp{- A,(6  +  Pi  -  Po)]  -  1 


and  by  a  similar  argument 


^#  {^)  =  [V’'(®.)1*M(1  “  P')(^  +  Pi  -  Po)  -  *  +  Po] 

(8) 

+  70Pt  +  7i(l  -  Pi)  +  o(l). 

where  p,-  denotes  the  right  hand  side  of  (7).  Substitution  of  (2),  (7),  and  (8)  into  (6)  proves 
(4);  and  (5)  may  be  proved  similarly. 


Remark.  If  we  compare  (for  example)  (2)  to  the  approximation  given  in  Theorem  10.16  of 
Siegmund  (1985),  the  new  approximation  is  often  only  marginally  better  numerically  than 
the  old  one.  Indeed,  in  the  case  of  a  normal  mean  —  7o  =  0  by  symmetry  and  the  two 
approximations  do  not  differ  numerically.  However,  the  error  o(l)  in  (2)  is  substantially  smaller 
than  the  o(b)  of  the  previous  result.  Consequently  it  seems  reasonable  to  conclude  that  the 
accuracy  observed  numerically  for  the  normal  case  in  Siegmund  (1985)  is  not  a  lucky  accident, 
but  can  be  expected  quite  generzdly.  For  the  exponential  distribution  as  discussed  in  Section  3, 
7i  -  7o  =  -17/18,  and  inclusion  of  this  term  converts  a  good  approximation  into  an  excellent 
one. 

To  obtain  an  approximation  for  E^(r  —  i/  |  r  >  »/)  we  first  proceed  heuristically  and  then 
more  rigorously.  Observe  that 

(9)  E^(r  -  1/  I  r  >  i/)  =  f  Ei  (t)P^{S^  -  minSk  G  dx  \  t  >  v}. 

J[o.h) 

According  to  Poliak  and  Siegmund  (1986),  as  v  and  6  — »  oo 

-  min5fc  Gdx\r  >  t/}  -*  P0„{maxSn  G  dx} 

k<t>  n>0 

in  the  sense  of  weak  convergence.  If  we  use  (4)  and  this  limiting  relation  in  (9),  we  obtain 

Eu{r  -  1/  I  r  >  I/)  “  [AiV’'(^i)r^{Ai(6  +  pi  -  Ef^Mi) 

(10) 

-  exp(AiPo)£^»o«*P(“^i^i)  +  +  “  Po)]}  +  7i. 

where  we  have  put  M\  =  maXn>o  Sn-  According  to  Siegmund  (1979,  Theorem  1) 

(11)  £^<o(A^i)  =  |Ao  ^1  -  Pi  + -|Ao|7i  +  o(Ao)  (Ao  — ►  0). 

By  a  similar  argument,  if  $q  and  9i  converge  to  0  at  the  same  rate 
E0O  exp(-AiMi)  =  (Ai  -  Ao)“Ml^o|  +  |Ao|Aipi 

-  “  7i)  +  o(Ai). 

Substitution  of  (11)  and  (12)  into  (10)  suggests  the  approximation,  for  >  0 
Ey(r  -  1/  \  r  >  i/)  St  [AiV'*(tfi)]“*{Ai(t  +  2pi  +  Aq  *) 

+  Aoexp[Ai(p]  +  po)]/(Ai  -  Ao)  +  exp[-Ai(6  +  Pi  -  po)]}  +  7i; 


(13) 


and  a  similar  calculation  with  =  0  yields  in  that  case 


(14)  E^{t  -  u  \  r  >  v)^  {b  + Pi-  po)*  -  2Ao’exp[Ao(pi  +  Po)]  +  71- 


In  Theorem  1  below  we  attempt  a  mathematically  rigorous  derivation  of  (13)  and  (14) .  We 
first  let  1/  and  6  — *  oo  holding  Oq  and  $i  fixed  and  obtain  an  approximation  for  Ev{t  -v  \  t  >  t/) 
which  contains  some  difficult  to  evaluate  constants.  We  then  let  $o  and  0i  converge  to  0  (at  the 
same  rate)  and  obtain  simple  expressions  for  these  constants  in  terms  of  pi,  po,  and  7i.  In  the 
case  =  0  this  calculation  leads  to  (14).  In  the  case  0i  >  0,  however,  the  second  exponential 
on  the  right  hand  side  of  (13)  does  not  survive  our  mathematically  rigorous  analysis,  because 
it  disappears  into  the  o(l)  error  term  when  6  — «  oo.  For  most  values  of  6  and  6i  the  term 
is  very  small,  suid  it  makes  little  difference  whether  it  is  included  in  the  final  approximation 
or  not;  but  for  small  values  of  $i  it  can  lead  to  a  significant  improvement.  See  the  numerical 
examples  of  Section  3. 


Let  Mo  =  minn>o  -SV,,  M  =  niaXn>o  Sn. 

Theorem  1.  Assume  Oq  <  0  <  Oi.  Suppose  i/  — ♦  oo  euid  6  -+  oo.  Then  for  $i>  0 
E^{r-i.\T>i^)  =  Eg, {S^^)/2Ee,iSr^)  -  Eg,{Mi) 

,  .  +  E|^{Mo)  [  <  oo}P(jo(Mi  €  dx) 

(13)  J[0,ool 

+  f  Ef  {Sr_\T-  <  oo)Pe„{Mi  €  dx)]  +  o(l), 

<110, oo) 

and  for  =  0 

Eu{r  -  1/  \  T  >  i/)  =  {b  +  Pi  -  Po)*  +  7i  ~  7o  -  p0oi^i)  +  [  €  dx) 

•llO.oo) 


(16) 


+  2po[Ei„{Mi)  -  f  ES{Sr.)P,,[Mi  e  dx)]  +  o(l). 

-llO.oo) 


If  in  (15),  $1  and  converge  to  0  in  a  fixed  ratio,  the  sum  of  the  four  terms  following  6  on  the 
right  hand  side  of  (15)  equals 

2pi  +  Aq^  +  Ao/[Ai(Ai  -  Ao)]  +  Ao(pi  +  po)/(Ai  -  Aq) 


(17) 


+  2^oAi(po  +  Pi)*/(Ai  -  Ao)  +  71^1  +  o(Ai). 


Similarly  if  $o  0  in  (16),  the  sum  of  the  terms  following  (b  +  pi  -  po)*  on  the  right  hand  side 


of  (16)  equals 


-2Ao^  -  2Aq^(po  +  Pi)  -  (po  +  Pi)*  +  71  +  o(l) 


Proof.  We  consider  only  (15)  and  (17).  The  proofs  of  (16)  and  (18)  we  similar.  To  prove  (15) 
it  sufRces  by  (9)  to  find  suitable  approximations  for  the  different  terms  on  the  right  hand  side 
of  (6)  and  integrate  these  with  respect  to  the  limiting  distribution  of  P„{5„  -  minfcc,/  St  €  dx  j 
r  >  i/}.  According  to  Poliak  and  Siegmund  (1986),  this  limiting  distribution  is  the  same  zis 
the  unconditional  limiting  distribution 

(19)  lim  P,/{St,  -  min  St  €  dx)  =  Pe.fMi  e  dx). 

u—oo  ^  0<k<v  '0  I  i  j 

By  the  definition  of  r,  Wald’s  identity,  and  (for  example)  Theorem  9.28  of  Siegmund  (1985), 
V'*(5i)P«i(r)  =  J5?«,(Sr)  =  E0^(St  -  minSt)  +  .B«,(minSt) 

(20) 

=  6  +  £;«j(5T+)/(2f;«,(ST^))  +  Ej,(Afo)  +  o(l). 

as  6  -+  oo.  Also  by  Wald’s  identity  and  the  renewal  theorem,  as  b  -*  oo  (recall  that  N  = 
inf{n:5n^[0,6)}) 

^'mEl(N)  =  {6  +  i?^(S*J/(2£;,,(5,J]}P,' >  6}  -  x 

(21) 

+  Ef^ {Sn\ Sif  <0)  +  o(1). 

Substituting  (20)  and  (21)  into  (6)  and  integrating  with  respect  to  (19)  gives  (15). 

Now  suppose  that  $o  and  converge  to  0  in  some  fixed  ratio.  It  follows  from  Theorem 
1  of  Siegmund  (1979)  and  its  proof  that 


^«i(^T+)/[2-E^»i(‘5r+)]  =  Pi  +  71^1  +  o(^i). 


(11)  holds,  and  similarly 


E$^{Mo)  =  -A,  ^  -  Po  -  ;:Ai7o  +  o(Ai). 


Hence  the  sum  of  the  four  terms  following  h  on  the  right  hand  side  of  (15)  equals 
2pi  +  Aq^  ~  (A7^  +  ^7oAi)  f  P$^{r-  <  oo}P«p(Mi  €  dx) 

(22)  >  j 

+  /  Ef  {Sr_  -  po;r_  <  oo)P(jo(Mi  6  dx)  +  7i(^i  +  rAo)  +  o(Ai). 


We  consider  the  first  integral  in  (22);  the  second  may  be  analyzed  similarly  and  is  slightly 
simpler.  Let  po{0)  =  E${S^_)/2Es{Sr_)  {0  <  0).  An  argument  along  the  lines  of  Theorem  2  of 
Siegmund  (1979),  but  with  somewhat  more  care  to  details  shows  that  for  some  K  >  Q 

(23)  {r_  <  co>  =  / 

I  exp[-Ai(i  -  po)][l  +  A,£;|.(S,_  -  Po(«0)  +  ^(A?)], 
where  o(Af)  and  0(Af)  are  uniform  in  x  and  0(e"^*)  is  uniform  in  Ai  provided  Ai  is  small 
enough.  If  we  split  the  range  of  integration  into  [0,  (logAJ"*)*)  and  [(log  A^*)*, +oo)  and  use 
the  first  approximation  in  (23)  for  large  x  and  the  second  for  small  i,  we  find  that 

-AJ"^  f  P/;(r-  <  oo)Ptfo(Mi  e  dx)  =  Et^[exp{-AiMi)] 

fOA\  •'[“•“I 


- 


‘  f  exp{-Aix)Ep[Sr.  -  Po{0*i)]P0o{Mi  e  dx)  +  o(Ai). 

•'(O.oo)  * 


The  method  of  proof  of  Theorem  1  of  Siegmund  (1979)  gives  (12),  and  by  a  similar  argument 

/  -  Po'iT-  <  oo)Pjj,(Mi  £  dx)  =  -Ai|Ao|7o/(Ai  -  Aq) 

(25)  ^ 

+  /  exp(-Aix)P^(5,_  -  Po{Ol)]P0o{Mi  6  dx)  +  o(Ai). 

’'[O.oo) 

It  is  easy  to  see  that  the  integrals  on  the  right  hand  sides  of  (24)  and  (25)  are  of  order  Ao,  and 
hence  when  (24)  and  (25)  are  added,  they  cancel  up  to  an  error  of  order  Aq.  (This  cancellation 
is  very  fortunate;  otherwise  the  final  answer  would  contain  some  constants  which  we  do  not 
know  how  to  evaluate.)  Substitution  of  (12),  (24),  and  (25)  into  (22)  and  simplification  yield 
(17). 


3.  Numerical  Examples. 

A  particul2Lrly  simple  and  important  special  case  of  the  approximations  of  Section  2  is 
the  normal  case:  dFo{x)  =  {2ir)~^f^ exp{-x^ /2)dx,  rl>{6)  =  6^/2,  Aj  =  20,,  and  as  pointed  out 


above 


Pi  =  -Po  —  -583,  7o  =  7i  =  7- 
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It  is  easy  to  compare  the  approximations  of  Proposition  1  with  values  determined  numerically 
by  Dobben  de  Bruyn  (1968)  or  Lucas  and  Crosier  (1982)  to  see  that  the  approximations  are 
very  accurate.  We  have  performed  a  Monte  Carlo  experiment  in  order  to  obtain  an  indication 


of  the  accur2u:y  of  the  approximations  in  (13)-(14).  The  results  are  given  in  Table  1.  The 
Monte  Carlo  experiment  had  2500  repetitions,  and  the  results  are  reported  ±  one  standard 
error.  Figures  given  in  parentheses  are  theoretical  calculations. 

The  theoretical  approximations  appear  to  be  excellent,  with  the  approximation  for  large 
v  being  very  good  even  for  v  =  10. 


Table  1 
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b  =  4.83,  jflol  =  .5  [E,,{r)  “  790) 
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The  approximations  of  Section  2  allow  us  to  make  the  following  simple  assessment  of  Lucas 
an  Crosier’s  (1982)  fast  initial  response  feature.  As  a  result  of  some  numerical  experimentation 
they  suggest  starting  the  CUSUM  process  at  i  =  6/2.  For  the  first  test  in  Table  1,  for 
6  =  4.87  and  x  =  2.44,  we  have  “  790.  Also,  for  9i  =  |tfo|  =  -S.  Eg^{r)  5.9  and 

E^{t  —  1/  \  t  >  u)  —  9.3  for  large  u.  Hence  we  have  obtained  about  a  40%  reduction  in  the 
average  delay  to  detect  an  immediate  change  at  the  price  of  a  1%  increase  in  the  average  delay 
to  detect  a  change  which  occurs  later.  It  follows  that  the  fast  initial  response  feature  with 
1  =  6/2  results  in  an  overall  savings  if  */  =  0  occurs  at  least  25%  of  the  time.  (This  conclusion 
changes  slightly  if  we  make  the  comparison  at  other  values  of  ffi.) 


A  second  class  of  important,  simple  examples  involves  detection  of  a  change  in  the  mean  of 


an  exponential  distribution,  or  what  is  almost  the  same,  a  change  in  the  intensity  of  a  Poisson 
process.  A  systematic  discussion  with  extensive  numerical  tables  is  given  by  Lucas  (1985). 
For  illustrative  purposes  consider  detection  of  a  change  of  the  parameter  A  of  an  exponential 
distribution  (A  =  mean~^)  from  Aq  to  AJ  >  Aq.  To  put  this  problem  in  the  canonical  form  of 
Section  2  we  set 

dFQ{x)  =exp(i  —  l)di,  x  <  I 

0  >  1, 

d  —  \j6  -  1,  and  V’(^)  =  6  -  log(l  +  0),  where  A  =  (AJ  -  Aq)/ log(A5/Ao).  In  terms  of  the 
exponentially  distributed  observations,  say  yi,y2,'  •  •,  we  have  n  =  1  -  Aj/t.  Note  that  ^  =  0 
corresponds  to  A  =  A. 

It  follows  easily  from  the  lack  of  memory  property  of  the  exponential  distribution  that 
the  Po-distribution  of  Sr_  is  unit  exponential  on  (-oo,0).  Hence  from  the  Wiener-Hopf  fac¬ 
torization  of  the  distributions  of  and  St_  (e.g.  Siegmund,  1985,  Theorem  8.41)  one  easily 
obtains  that  the  Po-distribution  of  Sr+  is  uniform  on  (0, 1).  It  follows  from  simple  calculations 
that  Pi  =  |,  7i  =  ^,  />o  =  —1,  and  70  =  1-  (If  instetwl  of  the  exponentially  distributed  y’s  one 
observes  continuously  a  Poisson  process  N{t),  so  S„  =  S(1  —  Ay*)  is  replaced  by  N{t)  —  At  in 
the  definition  of  r,  then  6  =  log(A/A)  and  by  continuity  po  =  7o  =  0-  If  our  test  is  to  detect  a 
decrease  in  A,  then  po  and  pi  are  interchanged,  m  are  70  and  71.) 

Table  2  gives  a  brief  comparison  of  the  approximations  of  Proposition  1  with  numeri¬ 
cally  determined  results  of  Lucas  (1985),  which  are  given  in  parentheses.  The  quality  of  the 
approximations  is  excellent. 

Recently  a  number  of  authors  (e.g.  Woodhall,  1983)  have  developed  algorithms  for  numer¬ 
ical  computation  of  the  entire  distribution  of  r.  Analytic  approximations  to  this  distribution 
will  be  discussed  in  a  future  paper. 


••q 
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Table  2 

Exponential  Distribution:  E^^t) 

(Ao  =  1,  A  =  1.25) 

Parenthetical  Entries  form  Lucas  (1985) 
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